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Abstract. The Uehling correction to the energy levels is presented in terms of the hypergeometric functions 
2-Fi. This presentation allows to derived various asymptotics and approximations. Further applications of 
this method to other atomic characteristics are also considered. 
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1 Introduction 

Vacuum polarization effects play an important role in quan 
turn electrodynamics of bound states, especially in the 
case of muonic and exotic atoms. They are responsible 
for the dominant QED correction to the Coulomb energy 
levels. Recent and planned activities on muonic hydro- 
gen [Illl] demand a precision theory, which in particular 
involves higher-order vacuum polarization contributions. 
For instance, for muonic hydrogen such contributions at 

■ the second and the third order were found in [US]. 

In this paper we consider contributions of the one- 
loop electronic vacuum polarization (the so-called Uehling 
potential). We intend to find such presentations of the 
Uehling correction, which could be successfully applied in 
higher-order calculations. 

The first-order vacuum-polarization correction to en- 
ergy levels in a hydrogen-like atom is the simplest QED 

' correction. With a known wave function it is usually not 

, a problem to perform a numerical computation for any 
desired state. However, there is a number of applications 
when an exact or approximate result in a closed analytic 
form is still preferable. Among them are the estimation of 
the correction for a number of levels (as a function of their 
quantum numbers), the calculation of a derivative of the 
energy shift (which is one of the approaches to find the re- 
lated correction to the g factor of a bound electron/muon 
[5]), estimations of relativistic and recoil corrections etc. 
The energy shift induced by the Uehling potential 

V v (r) = -^^[ 1 dv P ( V )e,J-^£=) , (1) 
7T r Jq \ V1-»V 
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where the dispersion density function is 

^(1-^/3) 
P( v ) = !_ v 2 > ( 2 ) 

is known in a closed form in a non-relativistic approxima- 
tion (see, e.g., [6]) as well as exactly in Za. The relativistic 
results for the Uehling correction to the energy were found 
for the circular states in an atom with the orbiting particle 
of an arbitrary mass m for spin 1/2 [7] and [8] exactly 
in Za. For an arbitrary state of a Dirac particle in the 
Coulomb potential the results were obtained in [9]. 

Here and throughout the paper the relativistic units in 
which h = c = 1 are applied; m e is the electron mass, Z 
is the nuclear charge and a is the fine structure constant. 

The analytic results expressed in terms of the basic 
integral Kb c (n n ) are cumbersome: 

Kbc{n n ) = K lbc (n n ) - -K 2 bc(K n ) , 
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where 3-F2 is the generalized hypergeometric function, 
B(x,y) is the beta-function and K n = Zam/(nm e ). 
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The generalized hypergeometric function 3F2 can be 
expanded as a series at low K n , but not at high values 
of K n . Even for k„ < 1 the hypergeometric series is not 
always appropriate because of the increase of coefficients 
for n>l. A real parameter of expansion for high n and 
low k„ is riK n rather than n n (see |10)h 

Due to all these problems it is hard to apply (J3]) di- 
rectly and, in fact, to find asympotics in [TllSj we used the 
integral presentation instead of the explicit expression: 



K bc (n n ) = dv 



p(v) 



KnVT 



(l_ w 2)6/2 



1 



(4) 



Here we describe another way for the calculation of 
K bc , which is free of all the mentioned problems with large 
n n and n and may be applied to the higher-order pertur- 
bation theory. In principle, the non-relativistic results for 
the Uehling potential can be expressed in terms of cer- 
tain elementary functions, which are still hard to use for 
asymptotics etc. The method developed here is equally 
efficient for non-relativistic and relativistic calculations. 



2 Calculation of upper and lower limits of the 
basic integral K bc (K n ) as model calculation 

To explain our approach, we simplify the problem and de- 
rive at first the upper and lower limits for K bc (K n ) instead 
of its calculation. 

First we change the variable in the integral presen- 
tation Q, introducing y = \/l — v 2 , and for the basic 
integral arrive to an equivalent form 



K bc (n n ) = K n / dy f(y)y c - b+1 ^1^(1 + Kn yy 
We note that the weight function 



(5) 



(0) 



monotonously increases from / m ; n = /(0) = 2/3 ~ 0.67 to 
/max = /(I) = V2 — 1.41. It is also fruitful to introduce 
a new kind of basic integrals: 



><- I d yy 

'0 



c-b+1 



(l- y y/2[l-(- Kn y)]-°. (7) 



Because of the monotonous behavior of the weight func- 
tion ^ with the rest of integrand in being a positive 
factor, we arrive to the upper and lower limits of K bc {K n ), 
by substituting f(y) in (|5|) by its minimal and maximal 
values, respectively: 



) < K bc ( Kn ) < V2Q bc ( 



(8) 



Furthermore, taking advantage of the simplification of the 
integral, we arrive in ([7]) to the generic integral presenta- 
tion of the hypergeometric function and find 



hc(Kn) =K C n B(c-b+ 2,3/2) 

x 2 Fi(c,c-b + 2,c-b + 7/2; 



We remind that in contrast to 3F2, the function 2F1 has 
many well-known properties and it is much easier to deal 
with. 

The ratio of the lower and upper limits in ([8]) is obvi- 
ously a constant \[2j3 ~ 0.5 for any values of b, c and n n . 
That is not a trivial issue by itself, because in the wide 
range from K n < 1 to K n > 1 the integral K bc {n n ) is 
changing by many orders of magnitude. We can get ad- 
vantage from the knowledge of the limits, e.g., by writing 
an estimation 



K bc ( Kn ) = (1.04 ±0.37) Q bc (K n ) 



(10) 



To take further advantages of the approach above we 
have three key problems to consider: 

— an improvement in the presentation for Q b c(K n ), which 
would be easier to apply for various k„; 

— a derivation of estimations similar to ([8]) and (fT0|) for 
an arbitrary state which involves a number of K bc in- 
tegrals with different values of b and c; 

— an improvement in the accuracy of the estimation (JTOj) . 

Below we proceed with all three problems and consider in 
addition further applications of the approach. 



3 Improved presentations for the basic 
integral Q bc ( 

K n) 

3.1 Low-/t„ adjusted presentation 

The presentation ([9]) for Q bc {n n ) allows to derive easily 
a series at low k„. Meanwhile, its applicability is rather 
limited since the coefficients b and c may be quite large 
in the case of n ^> 1 (see [10] for detail). As a first issue, 
we improve a presentation for Q bc (K n ) at low and medium 
K n . We achieve that by applying a proper transformation 
of the hypergeometric function, many of which are well 
known for (see, e.g., Appendix in [11), but not for 
3^2- In particular, applying Eq. (e,5) from [llj we find 



Qbc( K n) — 



1 + Kn 



B(c-b + 2,3/2) 



x 2 Fi c, 3/2, c- b + 7/2; 



1 + K n 



(11) 



(9) 



Apparently, the advantage of this presentation is that the 
argument of the hypergeometric function is always below 
unity (since n n is positive). We note the appearance of 
pre-factor (k„/(1 + K n )) c in (fTT|) - which carries most of 
the changes in the value of K bc (K n ), while the remaining 
factor is a slow-changing smooth function of K n (see Sect. 
6.3 in [TO]; cf. Eq. CGl in this paper). 

For an arbitrary state we discuss the problem in Sect.[4l 
and here we present examples for the non-relativistic case 
with the circular states (I = n — 1) only, for which the 
whole correction to the energy is determined by a single 
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Fig. 1. Comparison of the convergence of the hypergeometric 
function in relations ffi). and (|13[) for the Is state: (a) - the 
exact function Qbc(i^i), (b) - the first five terms of expansion 
of 2 Fi in 0, (c) - the same for (fTTj) . (d) - the same for (fT3|) . 



-ftThc integral [TlfTO] with 6 = 2 and c = 2rt. In this particu- 
lar case, we note another advantage of presentation (fTTj) . 
The convergency of the hypergeometric series 



2 Fi(a,a',f3;z) = 



r(a)r(a') 



E 

k=Q 



r(a + k)r(a' + k) z k 



r{i3 + k) 



is different in terms of K n for different presentations. While 
the series in (JSJ is finite only for a limited rangcQ of low 
K n , the series in the presentations (jTTJ) is convergent for 
< K n < oo. 

Indeed, the speed of the convergency of the series in 

and (fTT| is also quite different and in particular at low 
K n the convergency is much faster for high n for 2-F1 in 
(fTTj) than in ([9]). That is because in the former case for 
high n the parameters are: a ~ j3 ~ 2n and a' ~ 3/2, 
while in the latter presentation: a ~ a' ~ /3 ~ 2n. 

The results for Qbc( K n) related to the ground state 
and to a circular state at n = 10 are presented in Figs. Q] 
and [2 respectively. To check how fast the convergency is 
we consider only the first five terms of the hypergeometric 
series for each presentation of 2F1. 

The behavior of the five-term approximation in ([9]) is 
very notable (cf. also Fig. [7| . As explained in [TOjfTT] , the 
actual parameter of the naive 1ow-k„ expansion, based 
on ([9j), is nn n and at nn n — > 1 the error of the related 
five-term approximation becomes very high. In contrast 
to that, the presentation (JTTJ) allows an improved low-K n 
approximation which is successful at least up to tt n ~ 1. 

We can also check how many terms of the hypergeo- 
metric series we need to reach a proper level of accuracy. 
The results are collected in Figs. |3] and [4j One can note 
that the 5-term approximations are quite successful in a 
broad area. E.g., as seen in Figs.HHH t ne improved 1ow-k„ 
expansion is accurate up to K n ~ 1. 




Fig. 2. Comparison of the convergence of the hypergeometric 
function in relations ©, and (|13|) for the circular state 
with n — 10: (a) - the exact function Qbc(^io), (b) - the first 
five terms of expansion of 2-F1 in (O, (c) - the same for (|11[) . 
(d) - the same for (|13[1 . 





Fig. 3. Relative errors for the hypergeometric function expan- 
sions in relations (|ll|l (top) and (|13[) (bottom) as compared 
with the exact value for different number of non-trivial terms 
in the expansion (from bottom to up: 1,2,3,5 terms). The re- 
sults are for the Is state. 



1 E.g., for c > 1 and b 

CK n /2 < 1 [10]. 



2 the series is convergent for 
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3.2 High-K„ adjusted presentation 

In contrast to the presentation ([3]) for Kf, c with 3F2, the 
presentations (fTTj) for Qb c with 2F1 are easy to adjust for 
the high-K„ case. 

To derive a presentation efficient at high n n , we apply 
to (fTTj) the transformation (e,7) from [TT] and find 




Qbc{K"n) 



B(2-b, 3/2) 



x 2 F l c, 3/2, 6-1; 



1 + Kr. 



(12) 



{l + K n ) b - 2 B(-2 + b,c-b + 2) 



x 2 Fi 



b + 2, 7/2 -6, 3 -6; 



1 + Kr, 



We note that for the case of integer b, which actually 
takes place for any state, the expression contains a sin- 
gularity and should be properly regularized. In a general 
case we apply the substitute b — > b + e in (jlip , which is 
obviously valid, next apply the transformation (e,7) from 
[TT] and then consider a limit of e — > in l|12p . That de- 
livers us a finite well-determined result. Eventually, the 
result for b = 2 reads 



1 2c 



1 + « T . 



x <j 55-2*1 ( c - e ,--e,l-2e,— — 



ln(l + k„) - ^(3/2) - ^(c) + 2^(1) 



3 1 

X2.F1 I e,-,l,- 



2' ' l + K n 



(13) 



where VK^) is the logarithmic derivative of Euler's gamma 
function. A comparison of approximation of this presen- 
tation by a partial sum of the hypergeometric series is 
summarized in Figs. [JJ and [5J The results for general value 
of b is presented in App. [A] 

Obviously the high-/t„ asymptotics (cf. [7]) should con- 
tain In K n , which follows from well-known asymptotics of 
the vacuum polarization. The derived presentation allows 
to separate the logarithmic contribution explicitly. 

With appearing of the 2-F1 function of the argument 
1/(1 + K n ) one can consider the approximation of 2F1 by 
a few first terms of the hypergeometric series (see Figs. 
and[4|. We see that c = 2n and other parameters of 2F1 are 
of the order of unity, therefore for high n n the expansion 
is effectively done in n/n n (cf. [TO]). 

As long as we deal with © and (TTU|) , we can present the 
result for any state with combinations of Q2c{n n )- How- 
ever, improvements of estimation (|10|) . which we consider 
below, involve various cases of 2F1 and the transformation 
becomes more complicated (see App. lAlfor more detail). 




Fig. 4. Relative errors for the hypergeometric function expan- 
sions in relations ((TTJ (top) and (|13[) (bottom) as compared 
with the exact value for different number of non-trivial terms 
in the expansion (from bottom to up: 1,2,3,5 terms). The re- 
sults are for the circular state with n = 10. 



4 Estimations for corrections to the energy 
shifts 

The estimations above are achieved for the basic integral 
Kbc(^n), while the results for any state except circular 
ones contain a set of such integrals. In particular, the non- 
relativistic result for an arbitrary nl state is of the form 

ma 



AE nl = - 

TT 



a (Za) 2 m 



Fnl(K n ) 



F n l(K n ) 



(n + l)\ 



'-<i\(2l + i+l)\ 



(14) 



(15) 



2(ra r -i) 



where n r — n — I — 1 . 

This expression is a sum of terms with the same sign 
and we can apply to each Kb c separately the approach 
based on our approximation for f(y). Before doing that, 
we like to mention that the approach can be also applied 
to various other representation, such as a sign-alternating 
series [9] 



(16) 



i\(n r — i)\k\(n r — k)\ 

i=o fe=o y ' K ' 

(2/ + z + fc + l)! 
(2l + l + l)!(2l + k + iy K2 ' 2l+l+k+2{Kn} 



or a presentation with derivatives of Kbc [TOl 

n-l-l 



F n l(K n ) 



(n + l) 



1 



1 



n r !(2n-l)! ^ (2Z + i + l)! il 
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1 



2(n "- l) / n r \ V ( 2 _d_ 
\{n r -i)\) \ Kn dK n 

mi- i-u ( \ 2{l+l) F w {n n ) 



(17) 



which expresses the correction for an arbitrary non-rela- 
tivistic state in terms of the result for the Is state, which 
is well known (see two previous equations at n = 1, I = 0) 
[6]. One more presentation for F n i{n n ) can be found in 

Validity of estimations presented in Sect. [2] and simi- 
lar can be easily proved not through any explicit analytic 
presentations of F n i, but via related original expressions 
for the energy shift before any integrations are taken. 

Such an expression in both relativistic (Dirac's) and 
non-rclativistic case is of the form 



AE = / rf 3 



x P(r) - exp - _ 
r V VI 



r I dv p(v) 
o 

2Zam e r^ 



ir 



(18) 



Here v is the spectral parameter for the Schwinger presen- 
tation [T3] of the vacuum polarization, p(v) is the spectral 
function, which in the one-loop case is denned in ([2]), while 
P(r) is a positive weight function which describes the den- 
sity of the distribution of the electric charge of the atomic 
bound particle. The latter in the non-relativistic case reads 

Pnk(t) = W Ml 2 , 

where &(r) is the Schrodinger wave function of the related 
atomic state, while in the relativistic case it is denned as 

PnAr) = \f(r)\ 2 + \g(r)\ 2 , 

where /(r) and g(r) are the radial parts of the upper and 
lower components of the Dirac wave function. 

All evaluations in Sect. [2] are based on a certain ma- 
nipulation with p(v) and on the fact that the w-integrand 
is positive. The same can be seen in (fT8f . That means that 
for any state we can replace f(y) (which is a result of a 
certain transformation of p(v)) by its minimum and max- 
imum and in this way we arrive at the upper and lower 
limits for the whole F n i. To derive such limits we may 
apply any presentation for F n i listed above, by substitute 
Kb c in the right-hand side for its upper and lower limit. 

That would not be clear from the point of view of the 
presentations above by themselves, which include positive 
and negative contributions as well as derivatives. How- 
ever, substituting f(y) by /(0) (or /(l)), we should arrive 
at such a limitation. That is obvious, as seen from the con- 
sideration above for any presentation without derivatives 
(such as ([T5]) and (fT6|) ). Concerning (fTT]) . we remind that 
it is deduced from a presentation without derivatives by 
using certain recurrent relations [10] (see also [?])• Mean- 
while, the relations are maintained by the shape of j4| for 
any spectral function p(v) and thus allow any substitution 
for f(y). 



In the presentation (fl6|) b = 2 for all integrals and c can 
take various values, while in the case of (fT5|) the situation 
is opposite: c = 2n and b varies. The consideration above 
and in particular the high-K„ expansion in Sect. 13.21 is 
derived only for 6 = 2. Applying this approach to sum 
(fT6)) . one can arrive to more complicated identities (see 
App. [Al for more detail). 



5 Presentation of the weight function f(y) as 
infinite series exactly and its direct 
approximation by polynomials 

Above we found that simplifying the function f(y) in ([5]) 
(in particular, replacing it by a constant) we succeeded to 
simplify the integral and present the result in substantially 
simpler terms than in ([3]). 

In this section we consider a presentation of f(y) by 
infinite series in terms of either form 



f(y) = J2 Ckyk ' 

k=0 
oo 

/(y) = E c ^ 1 -y) fc - 



(19) 
(20) 



fe=0 



The finite series of this kind can already be used as an ap- 
proximation (see, e.g., Fig. [5j, but direct approximations 
lead to even simpler and more accurate presentations for 
the energy shift (see Sect. 15.21 below). 

Here we consider series and approximations for the ba- 
sic integral K} >c ( y n n ) 1 leaving discussions of uncertainties of 
partial sums of series for App. [B] 

Before considering any particular expansion or approx- 
imation based on IT9|) and ([20]) . we note that the weight 
function f{y) is positive as well as the integrand, as dis- 
cussed in Sect. 0J 

In particular, we start from 

AE = - [ dyf(y)R(y) , 



where R(y) is a certain positive function. Considering the 
approximation of f(y) by f(y) and introducing the cor- 
rection factor 

f(y) = N(y)f(y), 
we relate the approximation of the energy shift 



AE 



to the exact value as 



dyf(y)R(y) 



N(y) / dyf(y)R(y) = N(y)AE, (21) 



AE 



where y is a certain unknown intermediate value < y < 
1. To make a successful approximation it is enough to 
approximate the function f(y) in such a way that it does 
not exceed a certain margin: | iV(y) — 1 1 < 5. So, the relative 
uncertainty of the approximation is at least below S. That 
is valid for arbitrary R(y) and in particular for the energy 
shift for an arbitrary state in a muonic atom. 
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Fig. 5. Top: plots of f(y) ([5]) (bold line), its partial sums with 
first the first three terms for series (llQf) and ([20} (thin lines) 
and the linear approximation (|25|l (dashes). The maximal and 
minimal values of f(y) are denoted by the dotted lines. Bottom: 
the relative errors of the three-term partial sum (thin lines) and 
the linear approximation (dashes). 
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n 
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15 
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x k\ 


2<=+ 2 V^ r 

r(fe-5/2) 
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Table 1. Coefficients of the Taylor series for f(y) expansions 
(fj"9"l) and (|20]>. 



5.1 Presentation of /(y) by an infinite series and the 
results for the integral Kb c (K n ) 

Presenting f(y) as a Taylor series, one can find the coeffi- 
cients of either series (|19j) and (|20|) , which are summarized 
in Table [U 

Applying the series (|19[) to the integral ([5]) leads to an 
infinite sum 

oo 

K bc (K n ) = ^ Cfc Qb~k : c( K n) > (22) 
k=0 





fo 


h 


h 


h 


d 





-0.477333 


-0.418767 


-0.41490 


di 








-0.135733 


-0.15737 


d 2 











0.02412 


8 


13% 


0.7% 


0.03% 


0.003% 



Table 2. Coefficients of various approximations of f(y) in (|24[) , 
The approximations are denoted as fi{y)- 

where one can use for Qb-k,c & n appropriate presentation 
(see, e.g., ©, (HU), (USD etc.). 

A similar evaluation of ([20]) leads to the result 

x 2^1 (c,^ + k,c-b+ l + k; rr^—) -(23) 
V 2 2 1 + K n ) 

A presentation similar to © and (p~3|) can be also derived. 

An important property of the both series for f(y) is 
that the coefficients are quite regular. Their absolute val- 
ues decrease with k. While the coefficients in ([l~9"l) are reg- 
ularly sign-alternating (except of few first terms), coeffi- 
cients in (|20p are all negative (except of a few first terms). 
Such a regular structure allows a simple conservative es- 
timation of accuracy of a partial sum of either series (see 
App. [B]) . To conservatively estimate the remainder of the 
series, we do the estimation prior the y-integration, deal- 
ing only with a certain expansion of function f(y), which 
is indeed universal for any state. 



5.2 Polynomial approximations 

The series above allow in principle to reach any accuracy. 
However, one very seldom needs the accuracy substan- 
tially better than 0.01%, and it may be more fruitful not 
to expand the function, but to approximate it. 

To obtain an approximation for Kb c (n n ) we need to 
approximate successfully f(y), e.g., in a form 

f(y)^f(y) = f(0)-(l-y) + f(l)-y 

+ 2/(1 -j/)- (d + diy + d 2 y 2 + . . .) (24) 

and to tune the coefficients d k to minimize the value 

0<y<l 

For instance, a very rough linear approximation 

f(y) = My) = /(o) • (i - y) + /(i) • y (25) 

is already compatible with the partial three-term sums for 
either series considered above (see Fig. [5]). The results of 
various approximations are summarized in Table [21 

Including one more coefficient we improve the accu- 
racy substantially. Apparently, we can continue and reach 
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discuss possible further applications. Prior to the discus- 
sion we have to note that the application of some results, 
such as a presentation by the infinite series (|22|) and (|23l) , 
is quite straightforward. Some other need more considera- 
tion to estimate the fractional uncertainty. That concerns, 
in particular, the approximation (|26j) . When the complete 
y integrand does not change sign, we can apply an esti- 
mation from Table [21 otherwise we consider our approach 
as a 'good' approximation, the uncertainty of which is to 
be revisited in each particular case. 



6.1 Relativistic Uehling correction for a Dirac particle 




The relativistic expression for the energy shift of the re- 
state in a hydrogen-like atom with a Dirac particle (the 
electron, muon or antiproton) reads as [9] 



AE„ 



_a{Zaf 

nlj — o r nlj K^rt) 

7r n z 



(27) 



Fig. 6. Top: the relative error of the approximative functions 
^2,2n(ft„), corresponding to function f(y) = fi(y), for n = 1 
(solid line) and n = 100 (dots). Bottom: the same for function 

/(!/) = Ml/). 



any required accuracy. We already know that a success- 
ful polynomial approximation is possible because of the 
existence of series (fl9|) and (f20|) . meanwhile here we use 
a more direct strategy to minimize the uncertainty of an 
approximation with a fixed number of terms. 
Applying to ([5j> the approximation 



k=0 



where 



d = f(0), di = /(l)-/(0) + d Q) 
d k = tf fe _i - d k - 2 for k > 2 , 

we obtain the corresponding approximation for K^ c 

N 

K bc {n n ) = ^ dk Qb-kA K n) ■ (26) 

fc=0 

We expect our estimation of accuracy to be rather con- 
servative. A fractional error of the approximation f(y) is 
plotted in Fig.[S]for the non-relativistic basic integral K 2c , 
which, we remind, completely determines the energy shift 
for the circular state with n = c/2. 



6 Further applications 

Above we have developed a method for calculation of the 
non-relativistic Uehling correction to the energy. Below we 



F nlj (/^n) 



nV r(2C + < + !)<! 



where 



and 



(Za) 



X/ {](„.» 



Za 



(_!)*+* 



i,k=0 



il« -«)!*!«-*)! 
r(2C + i + k) 

r{2( + i + i)r{2( + k + i) 



x < m 



( v ) + (n — i)(n — k) 

\ 1) J 

i-k)\ 



1 - 



y - ( -i) i+ ' +1/2 (i + 1/2) 

V = - (Enij/m) 2 , 



(28) 



fx,, 



n - H , 
nr,n n /(Za) 



E n ij = m 



1 



{Zaf 



(C + <) 2 



-1/2 



is the exact relativistic energy of the state for the Dirac- 
Coulomb problem. The non-relativistic limit of the rela- 
tivistic expression above is given by (|16[) [9] . 

The efficiency of our approach in the relativistic case 
is illustrated for the 3d 5 / 2 state in Figs. [7]and[8l 

Our approach for the approximation consists of two 
parts: firstly, we approximate a combination of K bc by a 
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Fig. 7. Comparison of different approximations of (|29(l for 
3d 5 / 2 state: (a) - the relation (|29p . (b) - the first five terms of 
expansion for ((9j), (c) - the same for (d) - the same for 

(|13[) . Solid lines correspond to Z — 1, dotted lines to Z — 92. 



Fig. 8. Relative errors of different approximations of (|29[1 for 
3d 5 /2 state, (cl) - the first five terms of expansion for and 
K3 = 0.3, (c2) - the same for = 0.9, (d3) - the same for (|13l) 
and «3 = 8. The relative errors of the other approximations (o, 
b, c at K3 = 0.3, 0.9, 8, which are not shown, are above 40%. 



presentation with 2-F1 ■ 
(cf. (J251)) is 



The simplest of such presentations 



K bc {n n ) = /(0) Q bc {Kn) + [/(l) - /(0)] Qb-iA^n) ■ (29) 

Secondly, we approximate the latter by few first terms of 
the related 2-F1, which may have different arguments in 
different presentations and thus is well adjusted for ex- 
pansion in a certain region of k„. At both stages we can 
apply approximations of the same kind, varying the num- 
ber of approximation terms. That improves accuracy, but 
does not change situation in general, because we do not 
change the character of the expressions. 

As we have already proved in Sect.Hl all non-relativistic 
results for uncertainty due to the approximation at the 
first stage (see, e.g., Table [2]) are valid for the Dirac rela- 
tivistic consideration. 

To study the efficiency at the second step we rely on 
(|29p as an example. The functions Qbc{^n) are typical 
functions used in the approximations in Table [5] and we 
approximate Qbc, related to 3d 5 / 2 , by first five terms in 
HJ, (fTTj) and (fl3|) . The results are summarized in Figs. [7] 
and [5] and confirm that the efficiency of the approach does 
not depend on whether we proceed relativistically or not. 
Concerning the value of the relativistic effects in Figs. [7] 
and [SI we remind that the effects contribute to the pre- 
factors in ([28| . to the index c = i + k + 2£ of the Kb c (K n ) 



integral and to its argument k„. The figures present only 
approximations of K bc at a given argument k„, i.e. only a 
part of relativistic effects. 



6.2 Corrections to the bound g factor 

As found in [5j, to obtain the correction to the bound 
g factor in the case of an arbitrary potential for a Dirac or 
Schrodinger particle it is enough to know analytically the 
energy E at the related level of parametrical accuracy: 



ffbound(^j) 



2j(j + 1) 



l-2v 



ij 



di 



(30) 



Differentiating the function, for which we can numeri- 
cally control the accuracy of approximation, we cannot 
be sure that the numerical accuracy of the derivative is 
good enough. However, we expect that the better is the 
approximation of the energy the better is the result for 
the g factor. 



6.3 Relativistic Uehling correction for a Klein-Gordon 
particle 

Recently, a perturbative series of the Klein-Gordon bound 
particle was discussed [15] and the Uehling correction was 
found for circular states [8] . It was shown that the Uehling 
correction is still expressed in terms of Kb c (K' n ), where 



2n-2l- 1 
2n 2 (2l + l 



■{Za) 



(31) 



6.4 Non-relativistic corrections to the wave function at 
origin tf'(O) 

An accurate value of the non-relativistic wave function at 
origin, ^(0), was discussed in muonic and exotic atoms 
for a number of occasions. Rs value is important for the 
finitc-nuclear-size corrections in muonic and pionic atoms, 
for the hyperfine structure in muonic atoms [16] . for the 
pionium lifetime [TB1[T2] etc. 

The correction to the wave functions of the ns state 
^ns (0) induced by the Uehling potential is of the form 



5W ns (0) = / G; is (0,r) Vu(r) <Mr) d\ , 



(32) 



where G' ns {r' ,r) is the non-relativistic reduced Coulomb 
Green function. The values of G' ns (0, r) are known in a 
simple form (see, e.g., |17j ) and in particular 



Gi,(0,r) 



m e -Zamr , 

\ 2Zamr\ln(2Zamr) — ^(l)] 

2tt r I L J 

2(Z amr) 2 — hZamr — l| . (33) 
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The results are presented in terms well-defined at any 
arguments and using the well-known hypergeometric func- 
tion iF\. The approach allows to derive various approx- 
imations and find useful asymptotics. The free one-loop 
vacuum polarization is in principle not a problem for nu- 
merical calculations. Our concern is rather the higher or- 
der corrections and we expect our results can be applied 
there. 



10 10 10 10 10 10 1 




Fig. 9. The relative error of <5<?i s (0)/!?ii,(0) approximation as 
a function of m for / = fo (top) and f ~ fi (bottom) from 
Tabled 



For the Is wave function we obtain (cf. |16j 



gMg) 



n * Jo 



x (2 + + 3«iyi) + 2Kiy In 



1 + Kiy 



a 1 



:«i-K23(«i) 



■K 



2,2+e 



(34) 



Since the sign of the integrand in ([34]) does not change^], 
we use for its evaluation relations based on (|2ip . simi- 
lar to ones obtained above for the energy corrections. For 
instance, we can apply the simplest approximations for 
f(y) and apply conservative estimation of the uncertainty 
as presented in Table (see Fig. [5]). 



7 Summary 

Above we have developed a method which allows to study 
the Uehling correction to the energy levels and some other 
atomic characteristics. The method can be efficiently ap- 
plied to the non-relativistic consideration as well as in the 
relativistic case for Dirac and Klein-Gordon particles. 



We have found that it also takes place for the 2s state. 
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Appendix 

A Transformation of 2 Fi adjusted for the l/n n 
expansion 

The transformation (I12p contains singular terms and in 
the text of the paper we consider only the case of b = 2 
(see (fT3"l) ). For integer 6/2we apply the identity 

^ r(a + k)r{a> + k)r(/3) z k 



k=0 



r(a)r(a')r(p + k) k\ 



+ 



r(a + N)r(a' + N)r(p) z 



N 



r(a)r(a')r(0 + n) m 

x 3 F 2 (a + N,a' + N,l; (3 + N,N + l; z) . 

to separate the pole of the function for negative integer (3 
(we choose for that N = —/3+1). That makes the result for 
Qbci K n) more complicated, however, with z = 1/(1 + n n ) 
one can easily derive a related expansion for high n n . The 
result for integer b > 2 is 



1 



(-1) 6 



r(7/2-6) | 2(6-2)! 
1 



1 + «n 
2 F X ( c, 3/2, -1 + 6, 



1 

x [ln(l + n n ) - ^(3/2) - ip{c) + i>{b - 1) + ^(1)] 
d_ 
We 



F 2 [c-e, 3/2- e ,l; 



(1 + «„) 6 - 2 



r( c ) 



l-e,l 



b-3 

£r(2-6- 



e=0 



(6- 



fc=0 

fc-3)! ( I 



k)F\--b 



kl 



1 + «n 



(35) 



10 



E. Yu. Korzinin et al.: The Uehling correction in muonic atoms exactly in Za 



while for the case of integer 6 < 2 we obtain 

i \ (-i) b (i + Kn ) b ~ 2 r( c + 2~b) 



3(^71) 



1 + K n J r(c) 1 

a Fi [ c + 2 - b, 7/2 - 6, 3 - 6, 



(2-6)! 
1 



1 + K„ 

x [ln(l + Kn ) + ^(1) - ^(7/2 - 6) 
- V(c-& + 2) + V(3-6)] 

"de 



F 2 c + 2 -6 - e,7/2-6-e,l; 



3 — 6 — e, 1 — e; 



1 



E 



1 + Kr, 

(-i) fc r( c + fc)r(3/2 + fc) 

r(7/2 - 6) 



e=0- 



(1-6-fe)! / 1 



fe! 



1 + K r . 



(36) 



At 6 = 2 the summation in (1331) and 



vanishes 

(X)fc=o = Sfe=o = an< ^ both results coincide with (fT3"|) . 



B Estimation of uncertainty of series for /(y) 
and related results for the energy shifts 



Accordingly to (|21j) and the related discussion, to estimate 
the uncertainty of approximation of either a particular 
base integral or the correction to the energy for a certain 
level, it is enough to find accuracy of the approximation 
of f(y). Here we consider its approximation by a partial 
finite sum of series ([T9"f and (12"01) . 



B.l Finite y-series (|19[> and its uncertainty 

To estimate uncertainty of the presentation of f(y) by 
the sum of the first N terms of (|2"2"j) . we note that the 
coefficients of (fTSf satisfy the conditions 



|c fc+ i| < \c k \, 
c fe • c k+ i < 0, fc > 3 



(37) 



Therefore the fractional uncertainty of the sum is less than 
cjv//(1). E.g., for the first five terms of (j2"2"|) the relative 
error for Kb c is below 10%, and for the 14 terms it is below 
1%. 



and c' k < for fc > 3. That allows to estimate easily a sum 
of coefficients s' N = X^fcLjv c fc which enter the remaining 
series. The difference between the infinite series for the 
base integral and the TV-term partial sum has a fractional 
value which does not exceed s' N /f(l). E.g., the error of 
approximation of Kb c by the first three terms of ([23]) is 
below 20%, and by the first five terms the error is below 
1%. 
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B.2 Finite (1 — y)-series ( |20[> and its uncertainty 

To estimate uncertainty due to partial summation in (|23|) , 
we note that 

oo 

/M0) =£^4 (38) 

k=N 



